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ABSTRACT:  High-density  lipoprotein  (HDL)  plays  an  important  role  in  the 
transport  and  metabolism  of  cholesterol.  Mimics  of  HDL  are  being  explored  as 
potentially  powerful  therapeutic  agents  for  removing  excess  cholesterol  from  arterial 
plaques.  Gold  nanoparticles  (AuNPs)  functionalized  with  apolipoprotein  A-I  and  with 
the  lipids  l,2-dipalrnitoyl-sn-glycero-3-phosphocholine  and  1,2-dipalmitoyl-sn-glycero- 
3-phosphoethanolamine-N-[3-(2-pyridyldithio)propionate]  have  been  demonstrated 
to  be  robust  acceptors  of  cellular  cholesterol.  However,  detailed  structural  information 
about  this  functionalized  HDL  AuNP  is  still  lacking.  In  this  study,  we  have  used  X-ray 
photoelectron  spectroscopy  and  lecithin/ cholesterol  acyltransferase  activation  experi¬ 
ments  together  with  coarse-grained  and  all-atom  molecular  dynamics  simulations  to 
model  the  structure  and  cholesterol  uptake  properties  of  the  HDL  AuNP  construct.  By 
simulating  different  apolipoprotein-loaded  AuNPs,  we  find  that  lipids  are  oriented 
differently  in  regions  with  and  without  apoA-I.  We  also  show  that  in  this  functionalized 
HDL  AuNP,  the  distribution  of  cholesteryl  ester  maintains  a  reverse  concentration 
gradient  that  is  similar  to  the  gradient  found  in  native  HDL. 


KEYWORDS:  apolipoprotein  A-I ,  HDL,  gold  nanoparticle ,  cholesterol  cholesteryl  ester,  LCAT 


■  INTRODUCTION 

High-density  lipoprotein  (HDL)  is  involved  in  the  transport 
and  metabolism  of  phospholipids,  triglycerides,  and  cholesterol. 

A  mature  HDL  has  a  size  ranging  from  7—13  nm  in  diameter 
and  is  typically  composed  of  cholesteryl  ester,  triglycerides, 
phospholipids,  and  apolipoproteins.  HDL  is  also  known  as  the 
carrier  of  “good”  cholesterol  in  comparison  to  the  “bad” 
cholesterol  carrier,  low-density  lipoprotein,  as  it  is  believed  to 
remove  excess  cholesterol  from  arterial  plaques.  The  major 
apolipoprotein  component  of  HDL  is  apolipoprotein  A-I 
(apoA-l),  which  is  synthesized  in  the  liver  and  intestine  with 
a  mature  form  composed  of  243  residues.1  To  understand  the 
functionality  and  physiological  properties  of  apoA-I,  a  detailed 
knowledge  of  its  tertiary  structure  is  therefore  required.  The 
most  studied  apoA-I  structure  for  a  discoidal  HDL  is  a  “belt” 
structure:  lipids  are  encircled  by  two  antiparallel  apoA-I 
proteins.2-8  A  “trefoil”  model  for  a  spherical  HDL  has  also 
been  proposed  based  on  chemical  cross-linking  and  mass 
spectrometry  (CCL/MS)  studies.9  Because  of  the  dynamic 
nature  of  apoA-I  in  solution,  there  is  still  no  full-length  lipid- 
free  apoA-I  structure  available.  However,  several  structural 
models  of  lipid-free  apoA-I  have  been  proposed  based  on 
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CCL/MS  results  and  molecular  dynamics  (MD)  simula- 

..  10-13 

tions. 

With  the  importance  of  HDL,  numerous  academic 
laboratories  and  the  pharmaceutical  industry  have  tried  to 
develop  HDL  mimics  for  various  therapeutic  purposes.14,15 
Studies  have  demonstrated  that  functional  mimics  of  HDL  are 
able  to  bind  and  transport  cholesterol.16-18  Particularly, 
Thaxton  and  co-workers  have  employed  a  5— 6  nm  diameter 
gold  nanoparticle  (AuNP)  as  a  size  and  shape  controllable 
scaffold  to  assemble  l,2-dipalmitoyl-sn-glycero-3-phosphocho- 
line  (DPPC),  l,2-dipalmitoyl-sn-glycero-3-phospho- 
ethanolamine-N-[3-(2-pyridyldithio)propionate]  (PDP  PE), 
and  apoA-I  to  form  AuNP  templated  HDL-like  nanoparticles 
(HDL  AuNPs).16  A  follow-up  study  has  demonstrated  that  this 
lipid-functionalized  AuNP/protein  construct  is  a  robust  passive 
and  active  acceptor  of  cellular  cholesterol.1  HDL  AuNPs  have 
shown  potential  as  therapeutics  for  diseases  of  cholesterol 
overload,  lymphoma,  and  other  indications.17,19-22  AuNP 
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conjugates  are  known  to  have  high  tailorability  and 
biocompatibility,  which  makes  them  a  natural  candidate  for 
such  studies.23  Moreover,  AuNPs  can  be  detected  by  their 
optical  properties,  so  the  lipid-functionalized  AuNP  is  a 
promising  agent  for  monitoring  critical  human  performance 
parameters  of  HDL  such  as  cholesterol  transport,  cellular 
cholesterol  balance,  and  the  inflammatory  response.24,25 

To  further  optimize  the  functionalized  HDL  AuNP  for 
therapeutic  use,  it  is  essential  to  have  detailed  structural 
information  about  it.  Experimental  approaches  such  as  high 
resolution  transmission  electron  microscopy,  X-ray  photo¬ 
electron  spectroscopy  (XPS),  and  UV— vis  spectroscopy  have 
provided  valuable  insights  into  this  aspect,  yet  detailed  and 
localized  structural  information  concerning  the  construct  is  still 
lacking.26  Molecular  dynamics  simulation  provides  an  alter¬ 
native  approach  that  may  complement  the  current  experimental 
limitations.  In  this  study,  we  combined  coarse-grained  (CG) 
and  all-atom  (AA)  MD  simulations  to  simulate  the  self- 
assembly  of  lipids  and  cholesteryl  ester  on  the  AuNP/apoA-I 
construct  to  gain  insights  into  its  structural  properties,  and 
further  make  comparisons  with  XPS  and  functional  results. 
Because  the  self-assembly  of  lipids  on  AuNP  is  likely  to  have  a 
time  scale  longer  than  fis,  it  is  not  feasible  to  simulate  this 
process  directly  using  AA  MD.  Our  general  approach  here  is  to 
simulate  the  lipid  self-assembly  process  using  the  Martini  CG 
model,27-29  followed  by  switching  the  already  assembled  CG 
results  to  an  AA  model  and  performing  AA  simulations  for 
structural  analyses  that  can  be  compared  with  experiment.  0  We 
find  strong  agreement  between  theory  and  experiment,  and  our 
results  provide  detailed  insight  as  to  how  the  gold  HDL-like 
nanoparticle  construct  can  function  similarly  to  native  HDL  in 
its  uptake  of  cholesterol. 

■  RESULTS  AND  DISCUSSION 

Composition  of  HDL  AuNP  and  Simulation  Model. 

Experiments  have  previously  demonstrated  that  the  HDL 
AuNP  is  composed  of  a  5—6  nm  AuNP,  a  mixture  of  DPPC/ 
PDP  PE  lipids  — 175: 175),  and  ~3  apoA-I  proteins.16'17  It  is 
suggested  that  the  disulfide  bond  of  the  PDP  PE  lipid  is  broken 
during  adsorption  on  the  AuNP  surface,  which  results  in 
thiolates  of  MPDP  PE  and  2-MPT  (Scheme  Si).26  Because  of 
the  distinct  interactions  between  the  AuNP  surface  and  these 
lipids,  particularly  due  to  the  strong  Au— S  interaction  (~50 
kcal/mol),  it  is  expected  that  the  self-assembly  of  these  two 
types  of  lipid  on  the  AuNP  will  have  different  binding  patterns. 
Directly  simulating  this  HDL  AuNP  without  understanding  the 
pure  lipid  cases  may  lead  to  difficulty  in  interpreting  the 
resulting  morphology  as  there  are  many  factors  that  can 
contribute  to  the  final  structure.  Therefore,  we  started  with  a 
simple  system  consisting  of  only  pure  lipid  (plus  the  Au 
particle),  then  a  mixture  of  lipids,  and  finally  the  experimental 
HDL  AuNP  system.  We  placed  a  5  nm  AuNP  in  the  center  of  a 
210Ax  210Ax  210A  cubic  box  and  randomly  placed  the 
lipids  inside  this  cubic  box.  A  5  /is  CG  simulation  was 
performed  to  model  the  self-assembly  process.  After  the  lipids 
were  well  assembled  on  the  AuNP  surface,  we  backmapped  the 
CG  result  to  AA  model  and  performed  a  100  ns  AA  MD  for 
structural  analyses. 

Self-Assembly  of  Individual  DPPC  and  MPDP-PE/2- 
MPT  Lipids  on  AuNP.  We  first  attempted  to  simulate  the  self- 
assembly  process  at  300  K  with  350  pure  DPPC  lipids. 
Unexpectedly,  instead  of  binding  to  the  AuNP,  the  DPPC  lipids 
formed  a  bilayer  membrane  in  the  solution  (Figure  Si).  It  is 


possible  that  our  5  fis  simulation  time  is  not  enough  to  sample 
the  binding  process.  Another  possibility  is  that  the  DPPC  lipids 
are  in  a  gel  phase  at  300  K  that  makes  it  less  dynamic  to  diffuse 
to  the  AuNP  surface.  This  can  also  be  attributed  to  the  rigid 
AuNP  because  it  can  make  water  molecules  frozen  around  the 
solid  surface.  This  has  been  a  known  limitation  of  Martini  CG 
simulations.27  To  overcome  this  problem,  we  increased  the 
simulation  temperature.  When  the  simulation  temperature 
increased  to  323  K,  at  which  point  the  DPPC  lipids  are  in  a 
liquid  phase,  the  DPPC  lipids  were  able  to  self-assemble  on  the 
AuNP  within  5  /is.  They  wrapped  around  the  AuNP  and 
formed  a  donut-like  structure  with  a  diameter  ~140  A  (Figure 
1A).  The  lipid  tails  do  not  interdigitate  in  this  structure,  and  the 


Figure  1.  DPPC  self-assembled  on  AuNP.  (A)  350  DPPC  lipids 
formed  a  donut-like  conformation.  (B)  240  DPPC  lipids  formed  a 
trefoil  conformation.  AuNP  is  colored  in  yellow.  Headgroup  and  tail 
regions  are  colored  in  blue  and  cyan,  respectively. 


thickness  of  this  bilayer  is  ~50  A,  which  is  the  diameter  of  the 
AuNP.  We  also  performed  the  same  simulations  with  240  or 
160  DPPC  lipids.  Interestingly,  a  trefoil-like  structure  was 
found  in  the  case  of  240  DPPC  lipids  (Figure  IB).  When  the 
lipid  number  was  reduced  to  160,  a  donut-like  conformation 
formed  again,  with  a  similar  bilayer  thickness  but  a  smaller 
diameter  of  ~105  A  (Figure  S2).  In  the  donut-like  structure, 
the  DPPC  lipids  seemed  to  maintain  a  disc  bilayer  structure  if 
we  neglect  the  AuNP  region  (Figure  1A).  As  discoidal  HDL 
forms  a  similar  disc  bilayer  structure,2  this  suggests  that  AuNP 
can  serve  as  a  core  replacement  for  discoidal  HDL  and  the 
binding  of  the  apoA-I  on  the  AuNP  HDL  may  have  similar  belt 
or  trefoil  structure.  We  did  not  intensively  screen  different 
numbers  of  DPPC  lipids  in  this  study  since  our  primary  interest 
was  to  understand  systems  involving  lipid  mixtures  rather  than 
the  pure  DPPC  lipid  system.  As  a  result,  we  did  not  quantify 
the  ranges  in  DPPC  lipid  numbers  that  lead  to  donut-like  or 
trefoil-like  structures.  This  and  the  equivalent  studies  of  lipid 
mixtures  may  be  interesting  for  further  work. 

For  the  MPDP  PE/2-MPT  lipids,  because  of  the  strong 
interaction  between  the  sulfur  group  and  gold,  presumably  they 
will  self-assemble  preferentially  with  sulfur  coordinated  to  the 
AuNP  surface.  Indeed,  in  the  case  of  350  molecules  of  MPDP 
PE/2-MPT,  all  molecules  were  assembled  on  the  AuNP  within 
the  5  fis  simulation  time.  The  overall  morphology  was  like  a 
curved  bilayer  lipid  membrane  bound  on  the  AuNP  (Figure 
2A).  Most  of  the  lipids  in  the  inner  layer  were  laid  down  on  the 
AuNP  surface  and  loosely  packed  (Figure  2B),  while  most  of 
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Figure  2.  350  MPDP  PE/2-MPT  lipids  self-assembled  on  AuNP.  (A) 
Overall  morphology.  (B)  Inner  layer.  (C)  Outer  layer.  Headgroup  and 
tail  regions  are  colored  in  red  and  pink,  respectively.  2-MPT  is  colored 
in  purple. 


Figure  3.  Self-assembly  of  mixture  of  DPPC/MPDP  PE/2-MPT 
(175:175:175)  on  AuNP.  The  morphologies  of  lipids  on  AuNP  with 
(A)  all  three  types  of  lipids,  (B)  DPPC  lipid  only,  and  (C)  MPDP  PE 
lipid  only.  The  same  colors  as  in  Figures  1  and  2  are  used. 


the  lipids  in  the  outer  layer  were  oriented  perpendicularly  to 
the  AuNP  surface  (Figure  2C).  The  thicknesses  of  the  inner 
and  outer  layers  were  ~17  A  and  ~26  A,  respectively.  As 
expected,  the  2-MPT  molecules  were  located  mainly  on  the 
AuNP  surface.  In  addition,  no  belt  or  trefoil  conformations 
were  observed  among  different  lipid  numbers  for  MPDP  PE/2- 
MPT.  We  also  performed  the  same  simulation  with  PDP  PE 
lipid  (instead  of  MPDP  PE/2-MPT)  and  found  similar  layer 
thickness  and  morphologies  (Figure  S3). 

Self-Assembly  of  a  Mixture  of  DPPC  and  MPDP-PE/2- 
MPT  Lipids  on  AuNP.  Next,  we  simulated  a  mixture  of 
DPPC/MPDP  PE/2-MPT  lipids  (175:175:175)  on  the  AuNP. 
In  this  mixed  lipid  case,  the  MPDP  PE/2-MPT  molecules  seem 
to  dominate  the  overall  morphology.  The  two  types  of  lipids 
mix  together  and  form  a  bilayer  structure  on  the  AuNP  surface 
with  an  overall  morphology  similar  to  the  MPDP  PE/2-MPT 
result  described  above  (compare  Figures  3A  and  2A), 
suggesting  the  MPDP  PE/2-MPT  lipid  dominates  the  overall 
bilayer  structure  in  the  DPPC/MPDP  PE/2-MPT  mixture 
lipids.  This  is  attributed  to  the  sulfur  atoms  of  MPDP  PE/2- 
MPT  molecules  that  have  stronger  affinity  to  the  AuNP  surface. 
Moreover,  because  of  the  different  binding  strength  of  DPPC 
and  MPDP  PE/2-MPT  on  the  AuNP  surface,  it  is  expected  that 
their  mixture  will  form  an  asymmetric  bilayer  structure  on  the 
AuNP  with  the  inner  layer  having  the  higher  propensity  for  the 
MPDP  PE/2-MPT  lipid.  Indeed,  a  structural  analysis  showed 
that  ~66%  MPDP-PE  molecules  were  located  in  the  inner 
layer,  while  ~65%  DPPC  were  in  the  outer  layer.  Within  the 
same  lipid  layer,  these  two  types  of  lipids  were  randomly 
distributed  and  did  not  have  a  clear  clustering  pattern  (Figure 
3B,C). 

Self-Assembly  of  DPPC/MPDP  PE/2-MPT  on  Two  and 
Three  ApoA-l  Loaded  AuNP.  ApoA-I  protein  is  known  to 
have  a  well-defined  amphipathic  property,  providing  a  domain 
boundary  between  lipids  and  solution  that  helps  stabilize  the 
overall  structure  of  HDL.  The  experiment  demonstrated  that 
the  HDL  AuNP  construct  we  studied  contains  ~3  apoA-I.16 
Given  that  the  belt  and  trefoil  structures  of  apoA-I  on  HDL 
have  been  proposed  based  on  experiments,2-9  and  that  the 


HDL  AuNP  is  a  mimic  of  HDL,  we  assumed  that  in  the  cases  of 
two  and  three  apoA-I  in  the  HDL  AuNP  will  form  similar  belt 
and  trefoil  structures  on  AuNP,  respectively.  We  therefore 
preloaded  the  belt  or  trefoil  apoA-I  on  the  AuNP  construct 
(Figure  S4)  and  simulated  the  self-assembly  of  lipids  on  the 
AuNP. 

We  found  that  lipids  are  assembled  differently  on  the  regions 
with  or  without  apoA-I  protein.  Figure  4A  shows  the  overall 
morphology  of  the  three  apoA-I  loaded  case.  The  orientations 


Figure  4.  Self-assembly  of  a  mixture  of  DPPC/MPDP  PE/2-MPT 
(175:175:175)  on  three  apoA-I  loaded  AuNPs.  (A)  Overall 
morphology.  (B)  Type  I  orientation.  The  lipids  contact  apoA-I  with 
their  tails  and  lay  on  the  AuNP  surface.  (C)  Type  II  orientation.  The 
outer  layer  lipids  orient  perpendicularly  to  the  AuNP  surface,  while  the 
inner  layer  lipids  mostly  lay  on  the  AuNP  surface  and  pack  loosely. 
The  same  colors  as  in  Figures  1  and  2  are  used.  ApoA-I  is  colored  in 
gray. 
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of  lipids  have  two  distinct  types.  The  first  type  is  found  in  the 
region  under  apoA-I,  where  some  lipid  tails  contact  with  apoA-I 
or  the  AuNP  surface  and  the  rest  are  packed  in  between  them, 
which  results  in  a  horizontal  orientation  of  lipids  to  the  AuNP 
surface  (Figure  4B).  This  type  of  orientation  is  similar  to  the 
case  of  pure  DPPC  lipids  (Figure  1A).  The  second  type  is 
found  in  the  region  without  apoA-I  (Figure  4C),  and  the 
orientation  is  similar  to  the  mixture  lipids  case  (Figure  3).  In 
addition,  the  apoA-I  proteins  do  not  maintain  a  canonical  trefoil 
structure.  Instead,  the  trefoil  conformation  changes  to  a  T- 
shape  conformation,  with  a  large  sector  and  two  small  sectors. 
We  find  that  the  lipid  orientations  for  the  two  apoA-I  loaded 
results  are  similar  to  that  for  the  three  apoA-I  loaded  case 
(Figure  S5).  However,  the  two  apoA-I  loaded  case  have  more 
of  the  perpendicular  type  of  orientation  since  fewer  apoA-I  are 
in  the  HDL  AuNP. 

Distribution  of  Lipids  and  ApoA-I  on  AuNP.  Previous 
qualitative  XPS  data  suggested  that  the  outer  lipid  layer  of  HDL 
AuNP  is  composed  of  mixture  of  PDP  PE  and  DPPC  bound  to 
apoA-I.26  However,  the  detailed  composition  of  these  lipids 
within  outer  or  inner  layers  is  not  clear.  To  evaluate  the 
distribution  of  lipids  and  apoA-I  on  the  AuNP,  we  reanalyzed 
the  XPS  data  and  estimated  that  the  outer  lipid  layer  is 
composed  of  roughly  64—68%  DPPC  lipid  (see  Supporting 
Information  for  details).  To  make  comparison  to  the  XPS  data, 
we  calculated  the  radial  distribution  function  (RDF)  of  specific 
molecules  in  the  HDL  AuNP. 

As  mentioned  earlier,  the  molecules  that  contain  sulfur,  2- 
MPT  and  MPDP  PE,  prefer  to  stay  on  the  AuNP  surface. 
Indeed,  we  found  that  2-MPT  has  a  sharp  peak  at  ~28  A,  which 
is  in  the  AuNP  surface  region  (Figures  5  and  S6).  For  the 


Figure  5.  Radial  distribution  function  (RDF)  of  DPPC/MPDP  PE/2- 
MPT  (175:175:175)  on  three  apoA-I  loaded  HDL  AuNP. 

MPDP  PE  lipids,  the  peak  is  at  ~33  A,  which  also  indicates  that 
this  type  of  lipid  prefers  binding  to  the  AuNP  surface.  Since 
DPPC  lipids  have  weaker  interaction  with  the  AuNP,  it  is  not 
surprising  that  the  majority  of  them  are  not  located  on  the 
AuNP  surface.  By  taking  a  42  A  cutoff  (as  the  inner  layer  is  ~17 
A  thick)  to  separate  inner  and  outer  lipid  layers  for  the  above 
RDFs,  the  outer  layer  is  composed  of  67%  and  64%  DPPC  in 
the  two  and  three  apoA-I  loaded  AuNP  systems,  respectively.  In 
addition  to  the  above  RDF,  which  was  calculated  based  on 
whole  molecules,  we  also  calculated  RDFs  for  the  nitrogen 
atoms  of  DPPC  and  MPDP  PE  (Figure  S7),  which  provided  a 
more  quantitative  comparison  to  the  experimental  XPS  data. 
Table  1  shows  that  the  whole  molecule  and  nitrogen  atom 
RDFs  are  very  similar.  It  is  worth  noting  that  it  is  possible  that 
the  extraction  step  of  the  outer  layer  lipid  in  experiment  peels 
only  outer  layer  lipids  from  the  non-apoA-I  region.  In  other 
words,  lipids  in  the  space  between  the  AuNP  surface  and  apoA- 
I  are  protected  by  apoA-I  and  may  not  be  successfully  extracted. 


Table  1.  Percentage  of  DPPC  in  the  Outer  Layerc 


cholesteryl  ester 

DPPC  ( %)a 

dppc  (%y 

two  apoA-I 

0 

64 

67 

20 

67 

68 

100 

58 

63 

three  apoA-I 

0 

67 

74 

20 

63 

67 

100 

55 

U 

60 

aCalculated  from  whole-molecule  RDFs.  ^Calculated  from  N  atom 
RDFs.  CA  42  A  cutoff  in  RDFs  was  used  to  separate  inner  and  outer 
layer. 


In  that  case,  a  larger  discrepancy  between  the  simulation  and 
experiment  may  arise,  especially  in  high  apoA-I  loaded  AuNP. 

Distribution  of  Cholesteryl  Ester  on  the  Function¬ 
alized  HDL  AuNP.  Given  that  our  simulation  results  are  in  line 
with  the  experimental  lipid  distribution  for  the  AuNP  construct, 
we  next  considered  the  distribution  of  cholesteryl  ester  on  the 
construct  using  MD  simulations.  To  address  this  question,  we 
simulated  with  20  or  100  cholesteryl  esters,  which  account  for 
~1%  or  5%  volume  of  AuNP  HDL,  respectively.  We  simulated 
each  condition  with  either  two  and  three  apoA-I  using  the  same 
CG  and  AA  simulations,  mixing  cholesteryl  ester  molecules 
with  DPPC/MPDP  PE/2-MPT  (175:175:175)  lipids  from  the 
beginning  of  the  simulations. 

The  RDFs  from  constructs  that  consist  of  either  two  or  three 
apoA-I  loaded  with  lipids  and  20  cholesteryl  ester  molecules  on 
a  AuNP  revealed  that  the  cholesteryl  ester  molecules  were 
distributed  mainly  near  the  AuNP  surface,  which  was  near  the 
core  region  of  the  HDL  mimic  (Figures  6  and  S8).  For  the  rest 
of  the  molecules  in  the  construct,  their  distributions  were  very 
similar  to  what  we  found  without  cholesteryl  ester.  When  100 
cholesteryl  ester  molecules  were  loaded  onto  the  HDL  AuNP, 
they  were  still  mainly  found  near  the  core  region  of  the  AuNP 
construct,  however,  with  a  less  sharp  peak  in  the  RDF 
compared  to  when  only  20  molecules  were  included  (Figures 
S9  and  S10).  Past  experimental  work  has  suggested  that  the 
inner  core  of  a  native  HDL  has  a  higher  concentration  of 
cholesteryl  ester  and  that  the  packing  of  cholesteryl  ester  in  the 
particle  core  maintains  a  reverse  concentration  gradient  that 
can  assist  in  binding  more  free  cholesterol  to  the  surface  of 
HDL.1  Our  findings  on  cholesteryl  ester  binding  and  moving 
close  to  the  gold  nanoparticle  is  encouraging  and  suggests  that 
the  functionalized  HDL  AuNP  can  bind  cholesterol  in  a 
somewhat  similar  manner  to  native  HDL. 

Measurement  of  Cholesterol  and  Cholesteryl  Ester  in 
HDL  AuNPs  Incubated  with  LCAT.  Lecithin/cholesterol 
acyltransferase  (LCAT)  is  known  to  bind  HDL  particles,  is 
activated  by  apoA-I  and  converts  free  cholesterol  to  cholesteryl 
ester  that  is  then  sequestered  in  the  core  of  HDL.  As  modeling 
data  show  that  HDL  AuNPs  support  cholesteryl  ester  binding 
to  the  core  region  of  the  particle,  we  directly  measured  if  the 
HDL  AuNP  can  activate  LCAT  and  support  binding  of 
cholesteryl  ester.  First,  we  added  free  cholesterol  to  a  solution 
of  HDL  AuNPs.  Next,  we  added  purified  human  LCAT  into 
solution  and  incubated  the  mixture  overnight.  Using  the 
Amplex  Red  assay,  we  were  able  to  quantify  the  amount  of 
cholesterol  and  cholesteryl  ester  on  the  HDL  AuNPs.  By 
conducting  the  assay  with  and  without  cholesterol  esterase,  it  is 
possible  to  determine  the  contribution  of  cholesteryl  esters  to 
the  total  measured  amount  of  cholesterol  on  the  HDL  AuNPs. 
Using  this  method,  we  determined  that  the  HDL  AuNPs  have 
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Figure  6.  Self-assembly  of  a  mixture  of  DPPC/MPDP  PE/2-MPT/cholesteryl  ester  (175:175:175:20)  on  three  apoA-I  loaded  AuNPs.  (A)  Overall 
structure.  (B)  RDF.  Cholesteryl  ester  is  colored  in  blue. 
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Figure  7.  HDL  AuNP  support  the  activation  of  LCAT,  as  measured  through  Amplex  Red  based  fluorescence  measurement  of  cholesterol  and 
cholesteryl  ester.  (A)  In  the  presence  of  cholesterol  and  LCAT,  a  significant  amount  of  cholesteryl  ester  is  detected  (**,  p  <  0.005;  N  =  3  per  group). 
(B)  No  significant  amount  of  cholesteryl  ester  was  detected  using  HDL  AuNP  without  apoAI  (Bilayer  NP)  in  the  presence  of  cholesterol  and  LCAT, 
which  confirmed  the  specific  activity  of  LCAT  on  HDL  AuNPs  that  contain  apoA-I. 


both  free  and  esterified  cholesterol  after  incubation  with  LCAT 
(Figure  7A).  Further,  these  data  indicate  that  cholesteryl  ester 
makes  up  ~16%  of  the  total  cholesterol  bound  to  the  HDL 
AuNP.  As  a  control,  HDL  AuNP  without  apoA-I  (Bilayer  NP) 
incubated  with  cholesterol  and  LCAT  did  not  support 
cholesteryl  ester  formation  (Figure  7B),  which  demonstrated 
the  requirement  for  apoA-I  as  an  activator  of  LCAT  and  the 
formation  of  cholesteryl  ester.  The  modeling  data  support  that 
the  esterified  cholesterol  can  assume  a  more  central  location 
within  the  HDL  AuNP  particle. 

Simulation  Limitations  and  Future  Perspective.  In  this 
study,  we  have  assumed  that  the  apoA-I  forms  either  a  belt  or 
trefoil  structure  on  the  AuNP  as  these  two  structures  are  the 
most  commonly  studied  models.  However,  we  cannot  exclude 
the  possibility  that  other  apoA-I  structures  occur,  such  as 
double  superhelix, helix-dimer  hairpin,  or  horseshoe 
structures,34  and  it  is  not  possible  with  our  current  simulation 
technology  to  determine  which  one  is  correct.  Future  studies 
that  combine  both  theory  and  experiment  will  be  required  to 
further  address  this  question.  Another  limitation  is  that  since 
the  gold— sulfur  portion  of  the  force  field  is  a  simple  two-body 
potential,  it  might  not  be  able  to  fully  describe  the  influence  of 
facets,  edges  and  points  on  the  icosahedral  gold  particle.  In 
addition,  as  human  HDL  contains  different  compositions  of 
phospholipid,  the  study  of  different  combinations  and  types  of 
lipids  will  be  important  in  providing  critical  information  in 
diagnostic  and  therapeutic  applications  involving  HDL  AuNP. 


■  CONCLUSION 

In  this  study,  we  have  demonstrated  the  potential  of  computer 
simulation  in  modeling  a  previously  designed  AuNP-function- 
alized  HDL  construct  by  using  CG  and  AA  MD  simulations. 
On  the  basis  of  our  work,  and  as  a  starting  point,  we  found  that 
loading  160  and  350  pure  DPPC  lipids  on  the  5  nm  AuNP 
surface  led  to  a  donut-like  conformation,  while  a  trefoil-like 
conformation  was  found  in  the  case  of  240  DPPC  lipids.  This 
result  provides  comparison  to  natural  HDLs  where  the  apoA-I 
is  found  to  associate  in  discoidal  and  spherical  HDLs  in  similar 
fashion;  however,  this  system  is  not  as  directly  comparable  with 
the  HDL  AuNPs.  More  directly,  when  DPPC/MPDP  PE/2- 
MPT  (175:175:175)  were  loaded  on  the  AuNP  with  two  or 
three  apoA-I  proteins,  the  lipids  had  distinct  orientations  in  the 
regions  with  or  without  apoA-I.  In  the  region  with  apoA-I,  the 
lipids  orient  horizontally  to  the  AuNP  surface,  while  in  the 
region  without  apoA-I  the  outer  layer  lipids  orient  perpendic¬ 
ularly  to  the  AuNP  surface  and  the  inner  layer  lipids  pack 
loosely  or  lay  on  the  AuNP.  RDF  analyses  demonstrate  that 
molecules  containing  sulfur,  such  as  2-MPT  and  MPDP  PE, 
were  mainly  located  on  the  AuNP  surface.  In  addition,  the 
outer  layer  lipids  were  composed  of  64—74%  of  DPPC,  which 
is  in  close  correspondence  with  experimental  results  based  on 
quantitative  XPS  measurements.  Collectively,  these  data 
correspond  with  published  functional  data  demonstrating  that 
the  HDL  AuNP  supports  cholesterol  efflux  from  cells  by  way  of 
a  receptor  typically  accessed  by  lipid-poor  apoA-I,  namely  ATP- 
binding  cassette  transporter  A-l  (ABCAl),  and  also  through 
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receptors  more  associated  with  binding  lipid-rich  apoA-I 
associated  with  more  cholesterol-rich  HDLs,  like  ABCG1  and 
scavenger  receptor  type  B-l.  Finally,  we  showed  that  cholesteryl 
ester  molecules  are  mostly  in  close  proximity  to  the  AuNP 
surface  somewhat  analogous  to  the  core  region  of  native  HDL. 
The  cholesteryl  ester  distribution  and  functional  data  suggest 
that  the  HDL  AuNP  construct  is  able  to  maintain  a  surface-to- 
core  (low  to  high)  concentration  gradient  similar  to  the 
gradient  in  native  HDL.  Incubating  the  HDL  AuNP  with 
cholesterol  and  LCAT  demonstrate  the  HDL  AuNP  is  able  to 
efficiently  activate  LCAT  and  form  esterified  cholesterol  from 
free  cholesterol  and  phospholipid. 

Our  findings  have  demonstrated  the  value  and  great  potential 
of  computational  approaches  from  an  atomistic  model  in 
bridging  the  knowledge  gap  between  experimental  observations 
and,  particularly  in  this  case,  helping  us  to  better  describe  the 
functionalized  HDL  AuNP.  While  we  only  considered  belt  and 
trefoil  loaded  apoA-I  structures  in  this  study,  the  resulting 
simulation  protocols  can  be  applied  to  other  apoA-I  models  as 
well  as  to  different  lipid  types.  Further  study  of  different  lipid 
mixtures  as  well  as  different  apoA-I  models  will  be  interesting 
and  valuable  to  the  HDL  AuNP  field. 

■  METHODS 

Coarse-Grained  MD  Simulation.  On  the  basis  of  the  most  recent 
experimental  results,  the  shape  of  the  5  nm  AuNP  used  in  the 
experiments  is  icosahedral  rather  than  spherical  as  previously 
assumed.26  Moreover,  the  disulfide  bond  of  the  PDP  PE  lipid  is 
broken  during  adsorption,  which  results  in  thiolates  of  1,2-dipalmitoyl- 
src-glycero-3-phosphoethanolamine-N-3-mercaptopropionate  (denoted 
as  MPDP  PE)  and  2-mercaptopyridine  (denoted  as  2-MPT)  (Scheme 
Si).  Therefore,  MPDP  PE  and  2-MPT  molecules  were  used  in  the 
simulation  instead  of  the  PDP  PE  lipid.  The  Martini  v2.2  CG  model 
and  force  fields  were  used.27,28  MPDP  PE  and  2-MPT  force  fields  were 
parametrized  using  the  automartini  tools35  with  adjustments  to  match 
the  all  atom  MD  results.  Details  of  the  force  fields  are  listed  in  the 
Supporting  Information.  The  gold  force  field  was  taken  from  another 
study.36  A  5  nm  icosahedral  AuNP  was  built  using  a  script  from 
OpenMD37  software.  For  those  studies  having  apoA-I  present,  a  belt 
or  trefoil  structure  of  lipoprotein  apoA-I  was  placed  on  the  AuNP  to 
define  initial  structures  for  the  simulations.  Lipids  were  randomly 
placed  in  a  210  A  X  210Ax  210A  cubic  box  using  the  packmol 
package.  8  Each  system  was  solvated  with  ~7000  CG  water  molecules. 
10%  of  the  water  molecules  were  replaced  with  antifreeze  water 
molecules.  Na+  ions  were  added  to  neutralize  the  net  charge  of  the 
system.  All  CG  MD  simulations  were  carried  out  using  GROMACS 
4.5. 639  with  standard  simulation  parameters  associated  with  the 
Martini  force  field.  Simulations  were  performed  under  the  NPT 
ensemble.  Unless  specifically  mentioned,  the  simulation  temperature 
was  300  K.  A5//s  production  run  was  performed  for  each  case.  For 
the  apoA-I  loaded  cases,  to  prevent  the  apoA-I  from  collapsing  on  the 
AuNP  before  lipids  self-assembled  on  the  AuNP,  we  restrained  the 
apoA-I  structure  and  let  the  lipids  be  free  to  self-assemble  on  the 
AuNP  first.  This  was  done  by  applying  500  kj/mol  Cartesian  restraints 
to  the  protein  during  the  first  1  ps.  The  rest  of  the  production  run  was 
unrestrained  except  for  the  gold  atoms  which  were  always  fixed.  We 
performed  two  independent  runs  for  both  the  2  and  3  apoA-I  loaded 
DPPC/MPDP  PE/2-MPT / cholesteryl  ester  (175:175:175:100)  sys¬ 
tems,  the  most  complicated  systems  in  our  study,  and  found  that 
overall  morphologies  from  the  two  independent  runs  are  similar 
(Figure  Sll).  Therefore,  we  only  performed  one  run  for  the  remaining 
cases. 

All-Atom  MD  Simulation.  The  initial  structures  for  all-atom 
simulations  were  extracted  from  the  last  snapshot  of  the  5  fis  CG 
results.  The  backward  mapping  procedure  (CG  to  AA)  was  adapted 
from  other  studies.50  Each  system  was  neutralized  with  Na+  ions  and 
solvated  in  a  truncated  octahedral  box  with  ~56000  TIP3P  water 


molecules.  AMBER  ffl4SB,40  lipid  14, 41  and  GAFF42  force  fields  were 
used.  The  gold  force  field  was  adapted  from  earlier  work 43  AA  MD 
simulations  were  carried  out  using  AMBER  14.44  SHAKE  was  used  to 
constrain  bonds  to  hydrogen  45  The  particle  mesh  Ewald  method  was 
used  for  calculating  electrostatic  energy  with  a  10  A  nonbonded  cutoff. 
Each  system  was  equilibrated  as  follows.  First,  the  system  was 
minimized  for  5000  steps,  and  subsequently  heated  from  100  to  300  K 
in  500  ps  under  constant  volume  condition.  During  the  heating 
process,  10  kcal/mol/A2  Cartesian  restraints  were  applied  to  the 
proteins  and  lipids.  Next,  the  system  was  switched  to  a  NPT  ensemble 
at  1  atm  and  300  K  for  500  ps  with  the  same  restraints.  Then,  a  500  ps 
run  with  1  kcal/mol/A2  restraints  was  performed.  After  the  system  was 
equilibrated,  a  100  ns  unrestraint  production  run  was  performed  for 
each  system  under  the  NPT  condition  of  1  atm  and  300  K.  Gold 
atoms  were  restrained  during  the  whole  simulation.  To  test  the 
convergence,  we  measured  the  RMSD  of  the  apoA-I  structure  from  the 
3  apoA-I  loaded  DPPC/MPDP  PE/2-MPT/cholesteryl  ester 
(175:175:175:100)  system  and  found  that  protein  structure  is 
relatively  stable  after  50  ns  (Figure  S12).  Therefore,  the  last  50  ns 
results  were  taken  for  structural  analyses. 

Experimental  Characterization  of  HDL  AuNP.  The  XPS 
measurements  were  performed  as  described  in  an  earlier  paper.26  A 
detailed  analysis  of  intensity  results  was  performed  to  determine  the 
relative  populations  of  the  various  components.  Results  of  this  analysis 
are  given  in  the  Supporting  Information. 

Measurement  of  Cholesteryl  Ester  Formation  on  HDL  AuNP. 
The  HDL  AuNP  and  bilayer  NP  (HDL  AuNP  without  apoA-l)  were 
synthesized  as  previously  described.17  After  synthesis,  HDL  AuNPs 
and  bilayer  NPs  were  loaded  with  cholesterol  by  adding  10  //L  of  a  1 
mg/mL  stock  solution  of  cholesterol  in  ethanol  to  250  pL  of  100  nM 
NPs.  For  LCAT-treated  groups,  10  pL  of  a  10  nM  stock  solution  of 
LCAT  was  added  to  250  pL  of  NPs  and  allowed  to  incubate  overnight 
at  37  °C.  Following  overnight  incubation  of  NPs  with  cholesterol  and 
LCAT,  the  NPs  were  purified  by  centrifugation.  Free  cholesterol  was 
measured  using  an  Amplex  Red  Cholesterol  Assay  Kit  (Thermo  Fisher 
Scientific,  Waltham,  MA),  as  per  the  instructions  provided  by  the 
manufacturer.  For  cholesteryl  ester  measurements,  cholesterol  esterase 
was  omitted  from  the  test  reagents.  A  students  f-test  was  utilized  to 
analyze  data,  with  p  <  0.05  demonstrating  statistical  significance. 
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